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Electron transport through [normal metal]- [quantum dot]-[topological superconductor] junction 
is studied and reveals interlacing physics of Kondo correlations with two Majorana fermions bound 
state residing on the opposite edges of the topological superconductor. When the strength of the 
Majorana fermion coupling exceeds the temperature T, this combination of Kono-Majorana fermion 
physics can be observed: The usual peak of the temperature dependent zero biased conductance 
a(V — 0, T) splits and the conductance has a dip at T = 0. The height of the conductance side-peaks 
decreases with magnetic field. 

PACS numbers: 73.43.-f, 74.45.+C, 73.23.-b, 71.10. Pm. 
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Introduction: Recent theoretical investigation of topolog- 
ical materials reveals that topological superconductors 
can host Majorana fermions [l|48|- Specifically, Majo- 
rana bound states (MBS) reside at the ends of a 1-D 
topological superconductor (TS). Interest in the physics 
of Majorana quasiparticles is due to their non-Abelian 
statistics^. Hopefully, MBS can be realized on the 1- 
dimensional edge of a 2D quantum spin Hall insulator 
with proximity induced superconducting gap XI 



or at the ends of a one-dimensional semiconducting 
wire with spin-orbit coupling, in proximity with an s- 
wave superconductor 0,11]. In both cases a ID TS is 
formed. Two Majorana fermions can form a neutral 
Dirac fermion, but detecting it requires non-local mea- 
surements. Yet, it was suggested |12Hl4j that MBS can 
be probed in tunneling process, namely, a local mea- 
surement that is sensitive to interference between various 
MBS. 

Motivated by the above analysis, we consider electron 
tunneling through an N-QD-TS junction composed of 
normal metal lead, quantum dot and ID topological su- 
perconductor. As was already noted, the TS hosts two 
MBS on its ends [IKS]. Consequently, within a reason- 
able approximation, the tunneling problem is reduced to 
that of transport in an N-QD-MBS junction. Our interest 
is focused on the interrelation between the Kondo physics 
prevailing in the quantum dot and the MBS physics pre- 
vailing on the TS. The analysis is naturally divided into 
the weak (T > T K ) and strong (T < T K ) coupling 
regimes, where Tk is the Kondo temperature. The effect 
of a magnetic field acting on the dot (through a Zeeman 
term) is also analyzed since it is an important for ob- 
taining the ID TS. We start our analyzes with the weak 
coupling limit which can be studied by the perturbation 
theory 

The Hamiltonian of the junction 
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includes the following components: 1) Hq for the nor- 
mal metal lead, held at bias voltage V. 2) The Kondo 



part = + JTkm expressed in terms of dot- 

lead and dot-Majorana fermion exchange interactions 
Hlk and Hkm, obtained by applying the Schricffcr- 
Wolff transformation on the U/t — > oo Anderson Hamil- 
tonian with hopping energy t, level energy e and width 
T = 27r 2 |£| 2 AT(0) < |e|. 3) the Majorana term H M 
describes coupling of strength v between two Majorana 
fermions 1 on the left end and 2 on the right end of the 
TS . 4) The Zeeman energy BS Z of the dot subject to an 
external magnetic field B. 

Employing Nambu formalism in the 4 dimensional space 
[spin]® [electron- hole] , we have, assuming the dot is at 
x = 0, 



Hlk 
Hkm 



^ C t(0)Qr 2 c(0), 
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The definitions are as follows: ji (i = 1,2) are Ma- 
jorana fermion operators satisfying ji = jj, jf = 1. 
The coupling between two MBS is given by an antisym- 
metric 2 (g) 2 matrix v = v ( \ \ ). The Kondo cou- 
pling constant to the normal metal is Jt, = 2|i| 2 /|e|, 
and c(0) = {c^,c^,c^,—c\) T is the normal lead electron 
operator at x = 0. In Nambu representation, the cou- 
pling between the quantum dot and Majorana states is 
given by a vector % = (Af, Aj., AT, — A|) T . The Pauli r- 
matrices act on (c CT ; Aq-) and (cj.; A*) blocks. Below we 
put A-f = A, A^ = — iX [l6j|. Finally, the operator Q is the 
Numbu space extension of the exchange interaction 
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where s is the operator of electron spin. 
Keldysh actions and Green's functions (GF): The cur- 
rent operator is defined as the time derivative of num- 
ber operator of the normal metal lead J — edNr, /dt and 
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takes a form 



J =T^(0)QV 11 +H.C 
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Within the Keldysh perturbation technique, the normal 
metal free electron action reads, 




Slcad = n I dtYskCk 
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where g k 1 is the inverse GF for lead electrons. For later 
manipulations, the GF of the lead integrated over mo- 
mentum g = j-Sfe^fc is introduced, g 
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Here all entries are 4x4 diagonal matrices g (lu) = 
3 22 H = <JV(0)diag[(l - f(u - eV))(l, 1,0,0) + (1 - 
/(w + eV))(0, 0,1,1)] and g<H = -ff>(-w) = 
ii\T(0)diag[(/(w - eV)(l, 1, 0, 0) + f{u + eV))(0, 0, 1, 1)] 
where f(w) is the Fermi distribution function. We 
also use the GF matrix in rotated Keldysh basis g = 
Kgp z K~ Y (here Pauli matrices p act in Keldysh space 
and K = ( } - 1 ) /V2). 

To facilitate computation we consider the energy gap 
A of the superconductor as the highest energy scale in 
the problem. Then for low applied voltage eV < A in the 
weak coupling limit T 3> Tk only zero-energy Major ana 
operator projection of the total quasiparticle operator in 
the superconductor is important [l6l.fl8j. 

To average the product of dot spin operators we ex- 
press them in terms of mixed Dirac (/,/!) and Majo- 
rana (t] x , T) y , r\ z ) fermions [l7j|: S+ = rj z p; S- = fr] z ; 
S z = —ir) x r) y ; f = (r) x — ir] y )/\/2. In the diagrammatic 
representation of Fig. Q] below we use Keldysh GF for 
these fermions, 

1 



F\v) = 2mf{B)8(v — B), Ff{v) = 
' ' v — B + id 

1 



v + iS 
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The Majorana fermion action in Keldysh space follows 
from Eq.([2]) to zero order in tunneling to the dot is, 

S " = \ E / dtdt'Mt^Gl^tt')^'). (9) 



Here the inverse matrix GF , the Keldysh components of 
the matrix GF for Mojorana fermions at the ends of TS 
follow from Eq.([2]) are given by 



[G 



-[id t - 2i0], 



JmG n (u) = ImG H {uj) tanh — , 

G < (lu) = -2f(uj)iImG R (uj). 

Non-linear conductance: The calculation of the current 
diagrams (Fig. [T] ) is straightforward. The conductance 
a = 4j£ to the second order (Fig. [T]) upper row) is, 



(2) = aW(2u)(l + 3R(B)). 



(10) 






FIG. 1: Diagrams defining the contributions to the conductance 
due to Kondo interaction The upper panel represents the sec- 
ond order, while lower panel describes the third order contributions. 
Thick solid line - Majorana GF, thin solid lines - lead fermions GF, 
dashed and dotted lines - impurity spin. Dotted line is the Majo- 
rana F z GF, and dashed line is the Dirac fermion Ff GF. The left 
vertex in both sets has Keldysh index 1, other vertices have Pauli 
cr z matrix due to the doubling of the time contour in Keldysh tech- 
nique. 



where a = , and 

W(x) = I(cosh- 2 ^ 



cosh 



x — eV 
2T 



), (11) 



1. 2 B W(B + 2v) + W{B-2v) 

= 3 [C ° Sh 2T + ~ WW) + 

W(B + 2v) — W(B — 2v) tanh — ]. 

W(2u) 2T T l 

(12) 

In the Zero field limit (B -> 0) the function R(B) -> 1. 
There is a clear separation between potential and spin 
scattering. We also note that the temperature depen- 
dence of the non-linear conductance is quite distinct from 
the standard behavior of the quantum dot between nor- 
mal leads. This is due to the resonance form of the Ma- 
jorana GF. The third order contribution to the conduc- 
tance includes large logarithmic terms, that is the hall- 
mark of the Kondo effect in a tunneling system with topo- 
logical superconductor. Only these terms are retained, 



r(3) _ 



r 



= 3a[—]W(2u)K(B), 
ir\s\ 

where K(B) = k(B) + k(-B) and 

k(B) = l - C0S h- 2 ^L{2v + B) + 



(13) 



(L(2u) +L(2u-B))(l 



W{2v-B) 
W{2v) 

, B , v s 
tanh — tanh — , 

2T T' 
incorporates dominant logarithmic terms. 
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Here D denotes a high energy cut-off which corresponds 
to the band width in the normai iead. The two terms ( 
I10I13P are combined to yieid the total conductance in the 
weak coupling limit, 



a = aW(2v){l + 3[R(B) 



TK(B) 
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]} 



a e \lb) 



where a ex is given by the term in the square brackets 
multiples by 3. When B -> 0, K(B -> 0) = L(2u) and 
the total conductance becomes 

a = aWT2i/)[l+3(l + - T - r L(2i/))l. (16) 

7r|e| 

In Fig. [3] the total non-linear conductance is displayed 




FIG. 2: The total nonlinear conductance 1151 1 G = cr/a versus 
applied bias. The central high peak and lower plateau correspond 
to v = OAT, D = 1000T but different values of magnetic field: 
B=0 for the peak and B=3T for the plateau. The same for the 
dashed and dot dashed curves, the only difference is that v = 2.4T. 



For all four curves we assume 



=0.1. 



as function of the applied voltage. Two distinct peaks 
are resolved as the coupling energy v of the Majorana 
fermions exceeds the temperature, a fact which has nat- 
ural explanation. This is one of the central results of 
the present study since it combines the Kondo and Ma- 
jorana fcrmion physics, and encourages experimental ac- 
tivity in searching a realization of Majorana fermions. 
Under magnetic field the heights of the peaks decrease 
and reveals more complicated structure of cf(V). We can 
identify Kondo correlations also by studying the temper- 
ature dependence of the zero bias conductance G*. This 
is shown in Fig [3] We notice that G* in Fig [3] correlates 
with Fig.[2j for small temperatures where v > T the total 
zero biased conductance displays a dip, which is replaced 
by a sharp rise as temperature grows and exceeds v 
The strong coupling regime: At T < Tk we use the mean 
field slave boson approximation (MFSB) to estimate the 
zero bias tunneling conductance. As in the weak cou- 
pling limit we neglect the continuous spectrum above 
the superconductor gap. In addition to the conditions 
A T, eV, Tk we assume the inequality V s < T where 
is the tunneling width into the superconductor. In 




FIG. 3: The zero biased conductance as function of temperature 
at zero magnetic field. G = <r' 2 ' /(aT/T) is the second order con- 
tribution (dashed line), while G' stands for the total conductance 
in the same units as G (solid line). We take v = 0.4r and use value 
r/(27re)=0.1. 



this case contribution to conductance that originates 
from the continuous spectrum is proportional (at reso- 
nance) to 



r 



(17) 



that is very small indeed [l9|. The Kondo temperature 
Tk = D cxp[— 7r|e|/r] is determined solely by interaction 
with the normal metal lead. Indeed, in the weak coupling 
limit (T ^> Tk) MBS effectively represent the topological 
superconductor. The poor-man scaling renormalization 
group equations in this case read, 



dJ L 
d\nD 
<Mlr 
d\nD 



= -N(0)Jl 
= -N(0)J L J LR , 



(18) 



where Jlr = 2i|A|/|e|. The solution of these equations 
defines the above expression for the Kondo temperature 
and the conductance, 



G 



1 



peak 



A 2 W(2u + B) 
Uf N(0)T In 2 (fM) 



(19) 



where 



d(u,B) = ({T 2 + Av 2 ){T 2 + (2v + Bf 



1/4 



Actually the RG equations ([18]) show that Kondo insta- 
bility is related to the normal lead while the impact of res- 
onance tunneling through the MBS on Tk is irrelevant. 
We also notice that the scaling invariance of the exchange 
part of the conductance da ex / din D = is in agreement 
with the second RG equation ([HI). Thus if T < T K , then, 
unlike the case T ^> Tk, the condition A ^> Tk has to be 
supplemented by the inequality Ts < T in order to jus- 
tify the small contribution of the continuous spectrum of 
the superconductor. 
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To analyze the strong coupling limit of the N-QD-MBS 
system, we recall the Anderson hamiltonian for the quan- 
tum dot. In Nambu space the hamiltonian and the cur- 
rent operator acquire the form 

Htm = ^T z V 7l +tST z c(0)] + h.c. (20) 
H d = -[e($T z d + Un dt n di ], I = ^id f c(0) + h.c. 

In the limit U/t — > oo, employing the slave boson 
technique, the dot is empty or singly occupied. The 
elctron creation operator is written as, d\ = f^b where 
the slave fermion /J and the slave boson b mimic the 
singly occupied and empty dot states. They fulfill the 
constraint J2 a dlda + b^b = 1 that is encoded by in- 
cluding a Lagrange multiplier ?y the action S and replac- 
ing the U term in the dot Hamiltonian. At the mean 
field level the constraint is satisfied only on the aver- 
age. In the first step, a formula for the average current 
using the effective action S e ff is obtained. This action 
depends on two c- number parameters: the boson field 
b Cl and the chemical potential r\ c . In addition, the ac- 
tion is a function of one quantum source field 9 q which 
describes the interaction with the current /. Keldysh 
technique is employed to calculate the partition function 
Z(6q) = J D(Sddcrf)exp[iS] that depends on 6 q : Tak- 
ing the variation of In Z on 9 q — > we get expressions for 
the average current. I = (e/2h)5\nZ/56 q . After replac- 
ing the d>, d with (f\ f)bo and introducing the constrain 
into the action, the problem (in MFSB approximation) 
becomes gaussian. Integrating out the field of the normal 
lead variables, the quantum dot slave fermions /, f' and 
Majorana fermions 7 we are left with the fermion part of 
the partition function at 

\nZ f (6 q ) = trln^Gj 1 - O q Tb 2 T z [g, p x })}. (21) 

where GJ 1 = GJq - b^T is the inverse dot total GF. We 
keep only first order in quantum field 9 q which is enough 
to calculate current. Here Gj and T are matrices of the 
same form in Keldysh space: f = ( r " r r A ^ with entries 
as retarded, advanced and Keldysh components. The dot 
inverse GF GJq R = oq ® (id t — er z ) describes the nonin- 
teracting level with energy shift e —> e = e+rj c . The term 
r of the total quantum dot GF is the principal contribu- 
tion. We define T below and derive an expression for the 
average current. Differentiate InZ/ and performing the 
trace in Keldysh space we get a compact formula 

/ = Zl^ tr ^ G f T K G A + 2GflmT R Gfg K )T z }. (22) 

where ImT R = (T R — T A )/2i. The vertex T acquires a 
form 

f = Tg + Tz VG 11 p x V^T z . (23) 



Here G n = ( °^ G ^ J is the Majorana GF of 71 state. 
From Eq. ([23| we can obtain all entries 

2ImT R {uj) = -r + 2|A| 2 /mGf 1 (w)(l-fi) (24) 
T k {lu) = Tg K + 2i\X\ 2 (l - fi) tanh ^/mGfi 

where f2 = a y ® t z + a z ® t x + a x ® r y 
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FIG. 4: The zero biased conductance as a function of temperature 
at zero magnetic field in the strong coupling limit (T < Tk)- The 
dashed curve corresponds v = 0.5Tx and the solid line v = O.lTjf- 

The MFSB approximation is more reliable in equilib- 
rium where V — > 0. Therefore, we consider below the 
temperature dependence of the zero bias conductance for 
different couplings between two MBS. In equilibrium the 
mean field equations for boson field 60 and Lagrangian 
variable r\ q can be obtained by minimizing the free energy 
F on these two parameters. 

F = -TY^^HG^M+^bl (25) 

In Eq. (|23|) the last term is the slave boson part of the free 
anergy that is the result of constraint. The Matsubara 
GF of the dot G/ can be easily be obtained from G r (uj). 
Due to the condition T$ < T and Tk <C A the contribu- 
tion of continuous spectrum above the superconducting 
energy gap can be neglected. The mean field equations 
are 

b l = -T5^tr[T,G/(w n )] + l. (26) 
Vc = -T ^ tr[G/(w„)f (a; n )]. (27) 

The first equation |26|) fixes the level position e. In 
the Kondo regime the level is singly occupied e = 
or T) c = \e\. This approximate solution ij c of the sec- 
ond equation (|27p is used in the first equation to de- 
rive nontrivial solution (60 7^ 0) for the boson filed in 
terms of Kondo temperature. In the Kondo regime non- 
logarithmic terms in Eq([27|) can be discarded. Direct 
calculations show that such logarithms appear only for 
normal metal lead and the all terms related to the MBS 
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can be dropped. Thus we have b'^T/2 = Tk where Tk 
was defined above for the N lead in agreement with the 
weak coupling analyzis fig]) . For T < T K the Eq.(|22|) 
yields an expression for the linear conductance 



dx u 2 (x) + v 2 (x) 
G~o = 9 ~ T9 ~ X J cosh 2 x{x 2 +g 2 )A^Y 

2„,2/ 



(28) 



A(x) = x +{gr-g\v(x)) + g x u (x) - 2xg x v(x). 

where Go = 2e 2 /h, and the real functions u, v are 
defined from u(x) + iv(x) = x/[(x + iS) 2 — v 2 ]. Pertinent 
dimensionless parameters are: the interaction energy 
of Majorana fermions v = v/T, the tunneling rates 
g r = Tb 2 /AT = T K /2T, g x = 2\\\ 2 b 2 /T 2 . Figure 1 
displays the linear conductance versus temperature for 
two values the Majorana fcrmion coupling energy v and 
|A| 2 = TTk/8. The zero biased conductance as function 
of temperature displayed in Fig. [U reflects the similar 
dependence in the weak coupling limit T > Tk shown 
in Fig. [3] In both cases, if v < T, smearing starts at 
lower temperature and correspondingly the peak in the 
conductance shifts toward lower temperatures. 

Conclusion: Keldysh technique has been employed to cal- 
culate the linear and non-linear conductance in a system 
consisting of a quantum dot connected to a metal lead on 
one side and ID TS hosting Majorana bound states on 
the other side. The dot was tuned to the Kondo regime. 
Under certain approximations the whole system is re- 
duced to N-QD-MBS tunneling system in both the weak 
(T > Tk) and strong (T < T K ) coupling limits. The 
conductance has two peak structure if the coupling en- 
ergy v of MBS exceeds the temperature. Under a con- 
stant magnetic field, Ziman splitting occurs on the dot 
and reduces the heights of the peaks. The magnetic field 
may result in a more complicated structure of non-linear 
conductance peaks. Rcnormalization group analysis is 
performed in the weak coupling limit while the mean field 
slave boson approximation is used at T < Tk ■ Our anal- 
ysis shows that in an attempt to probe the features of 



MBS physics, the role of the Kondo effect is decisive. It 
is manifested by the occurrence of a strong temperature 
dependence of the zero bias conductance, and exposes 
distinct behavior of the non-linear conductance as com- 
pared with that of the simpler N-MBS tunnel junction. 
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